The Born-Oppenheimer Approach to the Matter- Gravity 

System and Unitarity 



C. Bertoni 

Istituto di Radioastronomia, via Gobetti 101, 4-0129 Bologna Italy 
F. Finelli, G. Venturi 

Dipartimento di Fisica, Universitd degli Studi di Bologna and I. N.F.N, 
via Irnerio, ^6 ~ 40126 Bologna - Italy 

February 7, 2008 



Abstract 

The Born-Oppenheimer approach to the matter-gravity system is illustrated in 
a simple minisuperspace model and the corrections to quantum field theory on a 
semiclassical background exhibited. Within such a context the unitary evolution for 
matter, in the absence of phenomena such as tunnelling or other instabilities, is ver- 
ified and compared with the results of other approaches. Lastly the simplifications 
associated with the use of adiabatic invariants to obtain the solution of the explicitly 
time dependent evolution equation for matter are evidenced. 







The Born-Oppenheimer (BO) approach has been extensively apphed to composite 
systems, such as molecules, which involve two mass, or time, scales [1]. Such an approach 
has also been suggested for the matter-gravity quantum system [2] in order to generalize 
the suggestion that matter follows semiclassical gravity adiabatically [3] (in the quantum 
mechanical sense). The plausibility of such an approach relies on the fact that the mass 
scale of gravity is the Planck mass which is much greater than that of normal matter. 
Thus one may consider the matter variables as the " fast" degrees of freedom whereas the 
"slow" ones are the gravitational variables. This is in analogy with the case of molecules 
where one considers the intermolecular distance to be the slow degree of freedom and the 
electron coordinates to be the fast ones. 

The purpose of this note is to briefly illustrate and compare the diverse approaches 
[4, 5, 6] with a particular emphasis on a possible violation of unitarity in the evolution of 
the matter system. For this it will be sufficient to consider a Friedmann-Robertson- Walker 
(FRW) line element: 

= _di2 + a^{t)gijdx'dx^ = a^{r]) (-djy^ + gijdx'dx^) (1) 

where Qij is the metric for a three-space of constant curvature k (which we shall always take 
equal to 1) and we have introduced a conformal time 77. We shall consider a matter-gravity 
action given by: 



dr/ 



']^m^a^ -rr?V + L^{a,cp) 



(2) 



where a dot denotes a derivative with respect to 77, m is the Planck mass, m^V the 
gravitational potential [7] and L*^ the matter {(p) Lagrangian which is allowed to depend 
on a but not a. Prom the above one obtains a classical Hamiltonian: 



H = -^+rr?V{a) + H'^{a,(t>) (3) 

where tTq = — m^d, is the matter Hamiltonian (which does not depend on tTq) and 
the classical Hamiltonian constraint is = 0. We note that the momentum constraint 
(diffeomorphism invariance on a space-like three-surface) is automatically satisfied in this 

minisuperspace model. 

On quantising one obtains a Wheeler-DeWitt (WDW) equation [8]: 
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(4) 



Sa^ 

where 'I' is a function of a and and describes both gravity and matter. Subsequentely 
one makes a BO decomposition of ^ as: 



*(a,0) = ^(a)x(a,(^) 



(5) 



where x(a, 0) is not further separable. Coupled equations of motions for x and ip may 
then be obtained by first substituing the above decomposition into eq. (4) and contracting 
with X* obtaining [2] : 



2m2 



{X\H^'\X) 

ixlx) 
1 



^^\{x\^ 

2m? (xlx) I da 



where we have introduced covariant derivatives: 



2m2 (xlx) 
{X\X)J da^^^r 



with [9]: 



^ d ^ - d 
T) = — + iA- V> = — ~iA 
oa oa 



(xlx) 



da 



and a scalar product: 



(xlx) = j d(f)x*{a,(l))x{a,<f>) 



where the integral is over the different matter modes. 

On now multiplying eq. (6) by x and subtracting it from eq. (4) one obtains: 



(6) 



(7) 



(8) 



(9) 



^ (^M _ ^^M^) ^ ^ |_ DX = - {&)) X 



^2 



d [ d d 



X 



(10) 
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and we note that the r.h.s. of eqs. (6) and (10) are related to fluctuations, that is they 
consist of an operator acting on a state minus its expectation value with respect to that 
state, and in the BO (or adiabatic) approximation they are neglected. 

In order to better understand our equations it is convenient to consider the following 
form for the gravitational wavefunction: 

^ ^ e-r^^»'v; ^ i-e-^r^'i-'Hs^'' (n) 

where N (which is related to the Van Vleek determinant [3, 4]) and S"^-^-' are real. On 
neglecting fluctuations in eq. (6), in the semiclassical limit for gravity S"^^^ will satisfy 
the following gravitational Hamilton- Jacobi equation [10]: 

which includes the back reaction of matter through the average matter Hamiltonian. We 
note = —m?a and in the classical limit also for matter one obtains the classical 

Einstein equation. 

On substituting into eq. (10) one obtains: 



drj J 



A, 



^ ±f{H«)dr,'-if''" da' A 



dlogN 
da 



B + {&)) 



X 



(13) 



where we have introduced the time derivative through ih-^ = —i^ ^^g^ '§a ■ 
on neglecting the r.h.s. of eq. (13), corresponding to the semiclassical limit for gravity 
and the adiabatic (or BO) approximation, one obtains the usual evolution equation for 
matter (Schwinger-Tomonaga or Schrodinger) . Thus it is natural to identify the matter 
wave function in the above as [11]: 



since in the above mentioned limits it becomes the usual Schrodinger wave function. 

One may now search for eventual violations of unitary evolution for matter by consid- 
ering: 
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X 



-i/"'"' da' (A- At) 



(xl^^^'lx) 



fi^ OlogiV 
m? da 



2m2 



X - c.c. > = 



(x|D|x) 



2m2 



(Xl (d^ - (D^)) Ix) 



c.c. ^ = 0, 



(15) 



where we observe that (D), ((D^ — (D^))) are zero and we have assumed (H^) is real. Thus 
unless one considers non Hermitian Hamiltonians (or the presence of tunneling phenomena 
leading to instabilities) there is no violation of unitarity. We note that in our approach, 
wherein we identified the gravitational and matter wave functions and equations of motion, 
one could envisage the presence of instabilities both in matter and gravitation and such 
that they compensate in the composite system [12]. 

Let us now compare our results with those of others. One approach [4] consists of 
expressing ^ as: 



e R 



(16) 



and expanding S in powers of m?: 



1 



S = m'So + Si + -^S2 + O 



9 "^-^ ' I A 



m^S'o(a) -\ Ka2{a) 



+ 



Si{a,(f)) H omia,' 



+ 



(17) 



This allows one to rewrite ^ in the following factorized form: 



\Nk 



(18) 



where Nk is also related to the Van Vleck determinant. We may now tentatively identify 
the above 'tpK and XK with our ^j; (eq. (11)) and x (eq- (14)) which satisfy: 



+ m^V + (ii™) V = -77—^ — 7^T^\ — V 



2m? da^ 



^rn? (xlx) 



(19) 



and: 
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X + 



2 da da 



m 



52 



I da"^ 



ixlx) 



X 



(20) 



Indeed, on substituing ipK and xk for -0 and x into eqs. (19) and (20) and retaining 
terms to 0{m'^), 0{m'^), 0{m~'^), one obtains from (19): 



(21) 



— hi 



Nk 2 



+ {H^)o = 



(22) 



N 



2 VnI 



(23) 



where a prime denotes a derivative witli respect to a. Eqs. (21), (22) and (23) are 



(H 



M\ 



(D2)o and {H 



respectively tlie equations 0{m^), 0{m^), 0{m^'^) and by 
we mean the corresponding terms of 0{mP) and 0{m^^) respectively. Analogously from 
eq. (20) one obtains: 



(24) 



to 0(771^) and 0(777, 2) respectively (where we have defined the c-number by H^xk = 
H^Xk)- On comparing the two 0{rnP) equations and eliminating {H^^)o one obtains: 

< + ^ - S'oS[ = (26) 

which agrees with the corresponding equation obtained on substituing for directly in 
the WDW equation (4) [4]. Similarly on composing the two 0{m''^) equations one can 
eliminate ((^^)_2 + ^(^2)o) obtaining: 

- SW2 + H%-S'oa', + i^h-^ = (27) 
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again in agreement with the corresponding expression obtained directly from the WDW 
equation [4]. Thus eqs. (26) and (27) confirm our identification of ipk and Xk with and 
X respectively. 

The above equations (24) and (25) then lead to the following expression for eq. (13): 



{( 



M 



drj 



XK 



XK + 



m^NK \Nk n 

f{x\&\x)o K A 



2m? 



ixlx) 



Nk 
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where we now have XKs = e ^ ^° r^^^ ^ ^^'^^ XK- Note that (x|D|x) 
= now becomes (see also eq. (24)): 



(28) 



{x\l\x) 



(xl^lx) 



d(^ 



XK 



0. 



(29) 



One may now check, as in eq. (15), whether a violation of unitarity occurs, obtaining 
on using eqs. (28) and (29): 



ih— 

OT] 



XKsXKsd(l) = iHSq J 



Xk^Xk 
oa 



c.c 



d(t) = 



(30) 



in agreement with our general result eq. (15). As before the result is a consequence of 
the presence of back-reaction terms {H^) and — — , indeed as can be seen from 
eq. (28) such terms are subtracted from the corresponding operators acting on xk, thus 
leading to zero on integrating over all matter configurations. The presence of back-reaction 
terms also modifies the definitions of Nk and a2 with respect to those employed elsewhere 
[4], indeed from eq. (22) and (23) one sees that these quantities are not associated with 
pure gravity, but include the mean effect of matter (see eq. (19)). 

One may alternatively attempt to use instead of Nk a prefactor Nq satisfying eq. (22) 
for {H^)o = (pure gravity) and introduce simultaneously an adiabatic phase factor in 
eq. (18), that is using instead of ipK- 



which will modify eq. (22) leading to [5, 13]: 



(31) 
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- mS', - (^)o) + ^ + {H^)o = inS',{^), + {H^), = (32) 

which however is in disagreement with the distinction between the adiabatic and dynamical 
phases [14, 16]. Indeed as a consequence of reparametrization invariance the dynamical 
phases of gravity and matter wave-functions cancel (Einstein equation). Similarly it is 
known that the light and heavy systems have equal and opposite adiabatic phases [1, 14]. 
Further the presence of the adiabatic phase and the fluctuations on the right hand sides 
of eqs. (6) and (10) are associated with corrections to the adiabatic approximation, that 
is they are related to transitions between different states and are expected to be small [15] 
(if such were not the case one would not have the usual evolution equations). 

Let us now compare our result with yet another approach [6]. Our equation (6) was 
obtained by substituing the BO decomposition into the WDW equation and contracting 
with X*. One may also expand |x) on a suitable orthonormal basis \l): 

|X) = E^K^)I0 (33) 
and contract with respect to (n|. One obtains [6] instead of eq. (6): 

^D^^ + mV+ ^ ' , o / I V V' (34) 



where: 



2w? " (n|x) 2m2 (n|x) 



d - d 

D„ = — + M„ ; D„ = — - lAn (35) 



with: 



("Ix) 

Of course if one multiplies eq. (34) by c* and sums over n eq. (6) is again obtained. 
As before, one may now multiply eq. (34) by x and subtract it from eq. (4) obtaining: 

V ('^Ix) J \ (n[x) J 

and on multiplying the above by c* and summing over n the usual result eq. (10) is 
obtained. Thus, again, we have separated the gravitational and matter equations of motion 
and equations (19) and (20) now become: 
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and: 



'fyM _ {n\H^\x) \ ^ T^d^dxn ^ _ {n\Bl\x) \ 

/ " 9a da 2rn? \d(^ (^^Ix) / '' 



where we have defined: 



The above expressions simphfy considerably if |x) = \vi) with Xn = (01'^)) in 
case one has: 



and: 



/"o-M / luMi \\ ~ , 91nV;„„9xn„ ( d'^ -2 I A . 

(iJ -(n|i/ l-)JXnn + ;^^^^ = -^(^^-(n|D„Jn)Jxn„ 

where we have correspondingly defined: 



with: 



_ d . - _ 9 . 

oa oa 



Ann = -i{n\^\n) = i{n\ — \n} 
oa da 



and: 



4'nn = e'r^--'"''A 



~ — „—i f Annda' 
Ann — An 



Similarly, instead of eq. (14) one obtains: 

where in this case dfiirj) is the trajectory obtained in the seniiclassical limit for ipnn- 

We now note that the r.h.s. of eq. (13) (similarly for eq. (43)) always consists of an 
operator acting on a state minus its expectation value with respect to that state, that is 
it is associated with fluctuations. Let us then denote by {v) (Xi^) the eigenstates obtained 
instead of |n) (xn) on omitting fluctuations in eq. (43). Correspondingly one will have a 
solution Xus to the Schrodinger equation: 



(i?^-i;i^)x.. = (49) 

Prom eq. (49) one can obtain some information about the orthonormal basis indeed 
on contracting the above equation with Xa ((^1) ^^^^ that it is identically satisfied for 

for (A I = {vl but for (A| / one has: 

(A| - ^n§^) = (50) 

The above equations (49) and (50) allow us to identify the |A) with the eigenstates of the 
time-dependent invariants [17] associated with our time dipendent matter hamiltonian 

On using eq. (14), in analogy with eq. (33), one may express Xs in our new basis: 

A 

which on substituing in eq. (13) leads to: 

V (-{H^) - MA + (A|F^|A) + hax{v)Axx - ih—) xx oc fluctuations (52) 
A ^ 

and on neglecting fluctuations (and denoting the quantities then obtained by the additional 
superscript 0) one has: 
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From the above one then obtains: 

X° = Y.cl{0)e-U'iMA^'W^.'-^r^a'A..^^ = E-a(0)xa. (54) 

A A 
which is the form for the general solution of the Schrodinger equation for a time dependent 

Hamiltonian in terms of the eigenfunctions of the time-dependent invariants [17]. 

If we do not neglect the fluctuations we may write: 

Xs = Ea [cl{0)+6,{r,)] e-U\m-m.'-^r^a'A,,^^ 

= J:x{cI{0) + 6x(.v))xxs (55) 

where the Sxirj) will be determined by the fluctuations and, on substituting into our general 
probability conservation constraint eq. (15), we obtain: 



A 

ihY: [cfmS.in) + Slir,)cl{0) + 5l{v)Sx{v) + Sl{v)Sxiri)] = (56) 

A 

which just reflects the fact that an increase of the weight S\{r]) of one state is compensated 
by the decrease of another. Thus there is no overall violation of unitarity, it is only on 
considering a single state that the evolution appears to be non unitary [16], and the result 
is obtained without having to resort to the limit — > oo[6]. 

Let us conclude: the approach we have illustrated consists of writing the total matter- 
gravity wave function in a factorized form ipx (Born-Oppenheimer) satisfying coupled 
equations, where ip describes gravitation with the back-reaction of matter and x describes 
matter on a curved background with both equations including fluctuations. It can be seen, 
on neglecting fluctuations and considering a semiclassical approximation for gravity, that 
the two wave functions ip and x can be interpreted as describing gravitation with a back- 
reaction due to the mean energy of matter and matter satisfying a Schwinger-Tomonaga 
equation and following gravitation adiabatically. The initial coupled equations themselves 
are a consequence of the factorization ansatz and do not involve any approximation. 
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Other approaches in general also write a factorized wave function which is substitued in 
the WDW equation, however they do not proceed to obtain coupled equations for matter 
and gravitation, but the diverse terms in their resulting single equation are identified by 
equating powers of rri^ and subsequentely a non-unitary evolution for matter is obtained. 
We have seen in one case [4] that if such an expansion is performed in our coupled equations 
which include the effects of all backreaction no difficulty arises with unitary evolution and 
in the other case [6] that there is no need to neglect lower powers of m? in order to obtain 
unitary evolution. 
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